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Global fits to all data of candidates for neutrino oscilla- 
tions are presented in the framework of a three-flavor model. 
The analysis excludes mass regions where the MSW effect is 
important for the solar neutrino problem. The best fit gives 
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and m 3 — m 2 ~ 1.11 eV indicating essentially maximal mix- 
ing between the two lightest neutrino mass eigenstates. 

PACS number(s): 14.60.Pq, 12.15.Hh, 14.60.Lm, 96.40.Tv 



I. INTRODUCTION 

Neutrino oscillations were reported last year by the 
Super-Kamiokande Collaboration Jl]] to have been con- 
vincingly seen by the Super-Kamiokande detector in 
Kamioka, Japan, in the study of the atmospheric neu- 
trino problem. This has rekindled interest to explain also 
other deviations from expected results in neutrino physics 
in terms of neutrino oscillations. There are at present es- 
sentially four different types of experiments looking for 
neutrino oscillations: solar neutrino experiments, atmo- 
spheric neutrino experiments, accelerator neutrino exper- 
iments, and reactor neutrino experiments. 

Most analyses of neutrino data so far have been per- 
formed in two-flavor models. Although this is quite il- 
lustrative, it can lead to wrong conclusions concerning 
the necessary number of degrees of freedom to explain 
data. It is, e.g., possible that the four neutrino scenario, 
advocated by several authors, is an artifact of this sim- 
plification. 

Recently, several papers ||-||] have been published per- 
taining at demonstrating the possibility to obtain an 
overall description in terms of a three-flavor model. Most 
of these investigations have focused on demonstrating the 
plausibility of the scenario, without truly fitting all the 
data. 

Here we would like to report on the first stage of a 
global analysis of all the neutrino oscillation data in a 
three-flavor model. We will assume that CP nonconser- 
vation is negligible at the present level of experimental ac- 
curacy. Thus, the Cabibbo-Kobayashi-Maskawa (CKM) 
mixing matrix for the neutrinos is real. There are then in 
principle five parameters that can be fitted in this model: 
Three mixing angles and two mass squared differences. 
(This is in contrast with two two-flavor models, which 
give only four parameters: Two mixing angles and two 
mass squared differences.) The two mass squared dif- 
ferences in the three-flavor model enter the argument of 



the sinodial function in combination with the ratio L/E, 
where L is the source-detector distance and E is the neu- 
trino energy. The argument of the sinodial functions re- 
lated to the energy dependence can then be treated as yet 
another angle. In this case, the probability function can 
be seen as given by a rotation of the flavor states by five 
angles. This is in analogy to the two-flavor case, which 
can be seen as a rotation by two angles, one of which is 
given by Am L/4E, Am 2 being the mass squared differ- 
ence of the two neutrino masses. 

In the present analysis, we will make the simplifying 
assumption that the relevant ranges of mass squared dif- 
ferences are such that the Mikheyev-Smirnov-Wolfenstein 
(MSW) effect in the Sun is negligible. This ex- 
cludes essentially the sensitivity region 10 -6 -10 -5 eV 2 . 
The absence of any day-night variation in the Super- 
Kamiokande and Kamiokande data is also consistent with 
neglecting the MSW resonance effect in the Earth for the 
electron-neutrino data. 

The differing results for the solar neutrino experiments 
present another problem. Since the Super-Kamiokande 
and Kamiokande experiments can resolve the energy for 
the solar neutrinos, these experiments are of a rather dif- 
ferent and more detailed kind than the radio-chemical 
experiments. It is possible to vary the disappearance 
rate by going to a mass squared realm, where the argu- 
ments, i.e., the angles that depend on the mass squared 
differences, are sensitive to the energy of the neutrinos, 
since this is different for the different detection meth- 
ods. However, in this case, the energy domains for the 
different experiments do not vary in the way expected, 
making the CI experiment the smallest (sa 0.33) and the 
Ga experiments the largest (« 0.6) of the expected rate. 
One would, in fact, expect the Super-Kamiokande exper- 
iment, with the highest mean energy for the neutrinos, 
to give the smallest result, whereas in reality it is in be- 
tween. 

Confronting this with the apparent lack of variation 
in the experimental solar neutrino detection probability 
with energy, at least up to approximately 13 MeV, in the 
Super-Kamiokande data, we will adopt the point of view 
that the mass squared differences are in a range that leads 
to the solar neutrino oscillations being averaged out to 
their classical values. By fitting the three mixing angles, 
we then account for the weighted mean of all the solar 
neutrino experiments. By nature of the sensitivity level 
of the accelerator and atmospheric neutrino experiments, 
we will choose the mass squared difference ranges as 
0.2 eV 2 < AM 2 < 2 eV 2 and 10~ 4 eV 2 < Am 2 < 10" 3 
eV 2 . Here AM 2 will regulate the accelerator experiments 
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and the multi-GeV atmospheric neutrino data, whereas 
Am 2 regulates the low energy atmospheric experiments. 
This choice is thus not inherited from two-flavor model 
analyses, but arises from the sensitivity of the different 
types of experiments. 

It is in principle possible that there would be a reso- 
nance effect for the atmospheric electron-neutrinos pass- 
ing through Earth. Our chosen ranges for the mass 
squared differences exclude this effect in the present fit, 
since Am 2 > 1CT 4 eV 2 . 

Since the reactor experiments are disappearance ex- 
periments with no signal seen, they can at best exclude 
certain ranges for the mass squared differences due to 
their sensitivity level. The two ranges of mass squared 
differences chosen here respect this sensitivity level. 

The paper is organized as follows. In Sec. n, we briefly 
discuss the formalism used in our analysis including nota- 



tions, Gaussian averaging, and minimization. In Sec. HI 



we discuss the choice of data from different experiments 
used in the minimization procedure. Sec. |y| describes the 
minimization procedure and includes the results thereof. 
We also show that our results are consistent with the ex- 
perimental data, assuming that all of the deviations from 
the expected values are indeed related to the physics of 
neutrino oscillations. Finally, Sec. |v| presents a discus- 
sion of the results and also our main conclusions. 



II. FORMALISM 



The probability of transition from a neutrino flavor 
state a to a neutrino flavor state (3 is given by 



a/3 



P(v a 
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where 5 a p is Kronccker's delta and Am 2 b 
The three mass squared differences Am^, Am^, and 
Am 2 3 are not linearly independent, since they satisfy the 
relation 



Amj! + Am 
Therefore, we will define 



Am 2 3 = 0. 
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AM 2 ee Am 2 2 = m 2 - m 2 , 
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and thus A 7 
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From Eq. (|4|) it is clear that P[pp — ► v a ) = Piv a — > 
vp) and since we have assumed that CP is conserved, i.e., 
the CKM mixing matrix U is real, it holds that P(vp — > 
i> a ) = P(y a — > vp). This implies that 



P{D a -» Vp) = P(l 



(8) 



A. Notations 



B. Gaussian averaging 



In the present analysis, we will use the plane-wave ap- 
proximation to describe neutrino oscillations. In this ap- 
proximation, a neutrino state \ v a ) with flavor a is a linear 
combination of neutrino mass eigenstates \v a ) such that 



v a) = ^2u* a \v a ), a = e,fi,T, 



(1) 



where the U a aS are entries in a unitary 3x3 matrix U. 
The unitary matrix U is given by 



Since in practice the neutrino wave is neither detected 
nor produced with sharp energy or with well defined 
propagation length, we have to average over the L/E 
dependence and other uncertainties in the detection and 
emission of the neutrino wave. We will here use the Gaus- 
sian average, which is defined by 



where 
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P(x)f(x) dx, 
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U = (Uaa) 



U^l C/ M 2 UfjS 
U T 1 U T 2 U T 3 



A convenient parametrization for U is 

C2C3 S3C2 s 2 

U= I -S3C1 - S1C3S2 C1C3-S1S2S3 sic 2 

S1S3 - S2C1C3 -S1C3 - S2S3C1 cic 2 



(2) 



(3) 



where St = sin^ and c% = cos^ for i — 1,2,3. This is 
the so called standard representation of the CKM mixing 
matrix. We have here put the CP phase equal to zero 
in the CKM matrix. This means that U* a = U aa for 
a = e, fi,r and a = 1,2,3. 
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Here I and 7 are the expectation value and standard de- 
viation, respectively. 

By taking the Gaussian average of Eq. (Q), we obtain 
the average transition probabilities from flavor v a to fla- 
vor vp as 

3 3 

(Pap) = S a p — 2 y^UaaUpgUabUpb 

a=l 6=1 
a<b 



1 - cos (2ZAm 2 6 ) e -27 2 (A m 2 ab ) 2 
a,f3 = e, /j, t. 



(10) 
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The physical interpretations of the parameters / and 7 
are the following: The parameter I deals with the sensi- 
tivity of an experiment and is given by I ~ 1.27 (L/E), 
where (L/E) should be given in km/GeV or m/MeV, 
and Am^ should be measured in eV 2 . Note that we 
will here use {L/E) = (L) / (E) . The parameter 7 is a so 
called damping factor. For large values of 7, the depen- 
dence on the mass squared differences will be completely 
washed out, since 1 — cos (2lAr 
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when 7 — * 00, and the transition probabilities {P a p) 
will just be dependent on the mixing angles 0j, where 
i = 1,2,3. This corresponds to the classical limit. In 
the other limit, 7 — > 0, we will just regain Eq. (Q) from 
Eq. (fig), i.e., lim 7 ^ (-Pa/3) = P a/3 . 



C. Minimization 

In order to obtain the mixing angles 9i, where i = 
1, 2, 3, and the mass squared differences Am 2 and AM 2 , 
we minimize the following object function 



F= F(9i,9 2 ,0 3 ,Am 2 ,AM 2 ) 

1 AW 
in 



pCKLp. 



(11) 



where w = X^i P w i ana - -^cxp is the number of neutrino 
oscillation experiments. The function F is a least square 
method error function for the relative errors of the neu- 
trino oscillation experiments with weights Wi > 0, where 
i = 1,2,... , -/V oxp . The quantities P; = (P a p)i and p? xp 
are the theoretical (model) and experimental transition 
probabilities for the ith experiment, respectively. The 
reason why we are using relative errors in the error func- 
tion and not absolute errors (as is usually done) is due 
to the fact that the different categories of experiments 
have rather different numerical ranges for the experimen- 
tal transition probabilities. In the rest of this paper, we 
will write P a /3, L, and E instead of (Pap), (L), and (E), 
respectively, to keep the notations simpler. 

In accordance with the discussion in the Introduction, 
the following constraints will be used in the minimization 
procedure 

< 0i < -, where i = 1,2,3, 

1CT 4 eV 2 < Am 2 < 1CT 3 eV 2 , 
0.2 eV 2 < AM 2 < 2 eV 2 . 



The weights in the object function ( |ll| ) will be chosen as 

1 

ne, ' 



(12) 



where n is the number of experiments used in the min- 
imization of the same category (solar, atmospheric, ac- 
celerator, or reactor) and u = APi/P° xp is the relative 
error in the ith experiment. 



Our choice for the weights is motivated by our wish 
to treat all categories of experiments on the same footing 
(the factor 1/n), and enhance experiments with small rel- 
ative errors (the factor 1/ci). Here again, relative errors 
are used, since the range of values of the probabilities 
are rather different. However, in all our cases, ti < 1. 
Higher powers of 1/ej would enhance the good experi- 
ments too much. Since AP; is a sum of the statistical 
and systematical errors, an underestimation of the latter 
will, of course, also overemphasize the relative weight of 
a certain experiment. 

For the damping factor of an experiment, we will use 



L I AL + AE\ 



E V L 



(13) 



where AL is a combination of the uncertainties coming 
from (a) the length of the detector and (b) the distance 
between the detector and the neutrino source and AP is 
the uncertainty in the mean neutrino energy. In Table |, 
we have listed all the experiments, which we have used 
in the minimization procedure. 



III. CHOICE OF DATA 

Most of the data used in our analysis are taken di- 
rectly from the quoted publications. (See Table | for the 
references to each experiment.) We have chosen in this 
first analysis not to discriminate between the various ex- 
periments in any other way than by their experimental 
uncertainties. 



A. Solar neutrino experiments 

For the chlorine (CI) and gallium (Ga) experiments, 
we will define the mean neutrino energy as 



(14) 



where S = . The index r labels the different solar 

neutrino sources. In the case of a CI experiment (Home- 
stake), r can be pep, 7 Be, 8 B, hep, 15 0, or 17 F, and sim- 
ilarly for a Ga experiment (SAGE or GALLEX), r can 
be pp, pep, 7 Be, 8 B, hep, 13 N, 15 0, or 17 F. The quanti- 
ties (E) r and S^ r ' are the average neutrino energy and 
standard solar model (SSM) prediction for the neutrino 
capture rates of the solar neutrino source r, respectively. 
Using the values given by the Bahcall-Basu-Pinsonneault 
1998 (BBP98) SSM HJgJI, we obtain the mean neutrino 
energy to be P Ga « 1.1 MeV for the SAGE and GALLEX 
experiments and Eci ~ 5.5 MeV for the Homestake ex- 
periment. 

As for the water-Cherenkov detector experiments 
(Super-Kamiokandc and Kamiokande), we have chosen 
the mean energy of the neutrinos to be P w _c ~ 10 MeV. 
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We will assume that the uncertainty in the energy is 
AE ~ E, respecting the threshold energy for detec- 
tion. For the Super-Kamiokande and Kamiokande ex- 
periments, the energy resolution is taken to be AE ~ 2 
McV. 

For solar neutrino experiments, the path length for the 
neutrinos is given by the mean Sun-Earth distance. The 
length of the detector, the deviation in the Earth's orbit 
around the Sun, and the solar radius are all small in this 
connection, giving AL/L ~ for solar neutrinos. 

B. Atmospheric neutrino experiments 

The path length of the atmospheric neutrinos is given 

by 

L = £,(<&) = \/R 2 cos 2 z9 + 2Rd + d 2 — R cos i9, (15) 

where •& is the so called zenith angle, R is the radius of 
the Earth, and d is the altitude of the production point of 
atmospheric neutrinos. We have taken R ps 6400 km and 
d ps 10 km, which gives L ps 1.22 x 10 4 km for the bin 
mean value cost? = —0.95 (MACRO and Kamiokande) 
and L ss 1.03 x 10 4 km for bin mean value cosi? = —0.8 
(Soudan, Super-Kamiokande, and 1MB). The bins — 1 < 
cos?? < -0.9 and -1 < cosi? < -0.6 have AL/L ~ 0.10 
and AL/L ~ 0.50, respectively. 

The atmospheric neutrino experiments present a spe- 
cial problem, since they are not entirely comparable to 
each other. The Super-Kamiokande experiment measures 
several parameters and has presented the "ratio of ra- 
tios" integrated over the energy (actually L/E) as does 
also the Kamiokande experiment. Other experiments can 
only measure the ratio of observed neutrino flux and ex- 
pected neutrino flux. The double ratio eliminates most 
of the uncertainty in the muon-neutrino flux. However, 
this function is also more complicated to fit, and has lost 
much of the information of path length and energy of 
the neutrinos. We have therefore used the directly mea- 
sured disappearance rates evaluated from the data for 
upward-goin g ev ents given in the tables and figures of 
Rcfs. |l] J7|,^4]-p7t . This makes it possible to use several of 
the other atmospheric neutrino experiments on the same 
footing. The disadvantage is the uncertainty in the ex- 
pected neutrino flux. This uncertainty is claimed to be 
of the order of 20% §. 

The main uncertainty of the neutrino flux comes from 
the uncertainty in the primary cosmic radiation flux giv- 
ing rise to the pions and kaons responsible for the muon 
and (later) electron neutrinos. Let the overall factor of 
flux uncertainty be denoted by rj a for neutrino flavor v a . 
Thus, e.g., the actual muon-neutrino flux is $ p = ?/ Al ( E > A1 , 
where $ M is the theoretical muon-neutrino flux calculated 
in Ref. |Q. Then, in the absence of neutrino oscilla- 
tions, the ratio of measured and theoretical flux is just 
Rp, = Vix- I n the presence of neutrino oscillations, on the 



other hand, we instead obtain 

_ <l e P eQ + I^P^ + $ T Pra n „v 

Ra = "T , (16) 

for any flavor a — e,fi, r. Here P a p is the probability 
of oscillation from flavor a to flavor (3. Since the ratio 
R = <f> A1 /<f> e is considered to be quite well determined 
theoretically, perhaps even to 5%, due to known branch- 
ing ratios of the pions and their decays, we conclude 
that the various factors rj a are essentially the same for 
a = e, \i. This is consistent with the picture, which says 
that the uncertainty is due to uncertainty in the primary 
cosmic ray flux. We will further also neglect contribu- 
tions from the r leptons in what follows. They should 
come from production of c or b mesons, but the relevant 
cross sections are small compared to the pion cross sec- 
tions. 

For the muon-neutrino ratio measured by several ex- 
periments, we obtain 

R, = V (P^ + ^Pe)j ■ (17) 

At E ps 10 GeV (Super-Kamiokande), we have R ~ 3 
PIP^I, and the last term can be neglected if P efl is 
small, which turns out to be the case. Thus, we obtain 
P^ ~ 77 _1 P M . We will therefore calculate the solution 
with the overall factor rj varying from 0.8 to 1.2 in the 
experimental values for R^ of upward-going events with 
cos i? < —0.6. 

For the electron- neutrino ratio R e — r/(P ee + RP^ e ) : 
the term RP^ e is not negligable compared to the term 
P ee , and can therefore not be neglected. Instead of using 
R e as input, we will analyze R e as an outcome of the op- 
timization process, i.e., in terms of the fitted parameters. 

Finally, the sensitivity for the atmospheric neutri- 
nos can be divided into two groups. The first group 
(Soudan and Super-Kamiokande) is sensitive to the small 
mass squared difference only, whereas the second group 
(MACRO, Kamiokande, and 1MB) is sensitive both to 
the large and (slightly) to the small mass squared dif- 
ferences. Changing the ratio from 0.8 to 1.2 times the 
quoted ones (Table |), will therefore influence the deter- 
mination of the mass squared differences. 

C. Accelerator neutrino experiments 

The mean energies and path lengths for the accelera- 
tor neutrino experiments are all given in Refs. p9|-|22]|. In 
most cases, a measure of AL is given by half the length 
of the detector, giving AL/L in the range 0.15-0.35. For 
the transition probabilities, we have used the presented 
90% C.L. upper bounds, except for the Liquid Scintilla- 
tor Neutrino Detector (LSND) experiment, for which we 
have used their quoted — > v e transition probability. 



4 



D. Reactor neutrino experiments 



for the mean value and 



The reactor neutrino experiments are all disappear- 
ance experiments with no signal detected. We have used 
the published 90% C.L. bounds on disappearance signals 
mainly to limit the mass squared ranges. Once these are 
outside the sensitivities of the reactor experiments, the 
solution will automatically reproduce these. For the re- 
actor experiments, the uncertainties in L are all small 
compared to the uncertainties in energy determination. 



IV. NUMERICAL ANALYSIS 

To be able to minimize Eq. ( |TT| ) we have to use nu- 
merical methods, since there is no successful way to do 
this with analytical methods, the present function be- 
ing to large and we have far too many variables. Our 
procedure will be as follows. First, we will use a very 
simple stochastic procedure. We choose some number 
(n) of random points in the specified domain and evalu- 
ate the object function (|li"|). We simply take the point 
which gives the smallest value of this function. To obtain 
good statistics we repeat this procedure N times. Includ- 
ing all 16 experiments, discussed in the previous section, 
and choosing n — 10 6 and N = 20, we obtained 

0! = 28.7° ±0.7°, 

2 = 4.4° ±0.1°, 

6» 3 = 45.1° ±0.3°, 

Am 2 = (2.89 ± 0.21) x 10~ 4 eV 2 , 

AM 2 = 1.12 ±0.04 eV 2 , 

where we have used 

1 

x — x ± — =o"jv— i. 
VN 

The quantities x and <jn-i are the sample mean value 
and sample standard deviation, respectively. The point 
which generated the smallest value of the object func- 
tion rtllf) in our simulation was 



Am 2 



9i«28.9°, 6> 2 ~4.2°, 6 3 « 45.0°, 
2.87 x 10" 4 eV 2 , AM 2 ps 1.11 eV 2 



01 w 27.97°, i 
Am 2 w 2.84 x 10~ 4 eV 2 , 



4.45 c 



AM 



re 43.92° 
1.14 eV 2 



Having these 20 data points, we can now use them as 
initial points in a deterministic minimization procedure. 
We will use a sequential quadratic programming (SQP) 
method. The results obtained are 



h = 28.7° ±0.5° 
9 2 = 4.2° ±0.1°, 
h = 45.0° ± 0.1 c 

,2 



Am 2 = (2.87 ±0.22) x 10 
AM 2 = 1.12 ±0.02 eV 2 



" 4 eV 2 , 



for the best point. 

Finally, we have evaluated the theoretically expected 
probabilities for the different experiments, using the best 
point solution, to demonstrate that they are indeed con- 
sistent with experimental data. This is shown in Table ||. 

Upon varying 77, we have made an analysis to see the 
effect on the different parameters in the model. The pa- 
rameters then vary as is shown in Fig. ^. The stability 
of 63 = 45° is to be noted. The mixing angle 82 varies 
only slowly around 5°. The mixing angle 8\ varies from 
15° up to 45° in the interval of r\ investigated. Just as 
expected, the mass squared differences change slowly to 
compensate for this change in the atmospheric values. As 
mentioned earlier, the reason for this is that the two first 
of the atmospheric experiments are sensitive to the small 
mass squared difference and the three others to the large 
mass squared difference. 



V. DISCUSSION AND CONCLUSIONS 

From the results in Table ||, we see that the parameter 
set given as best point is indeed consistent with most 
of the experiments. It is clear from the beginning that 
the solar neutrino experiments cannot all be reproduced 
with the parameter range used here, as was mentioned in 
the Introduction. The more detailed Super-Kamiokande 
and Kamiokande results, that include energy resolution 
of the solar neutrinos, are very well fitted by this solution. 
The GALLEX and Homestake experiments are not so 
well fitted. In the case of GALLEX, the fit is within 
two standard deviations, whereas the Homestake result 
is several standard deviations off. 

For the atmospheric neutrino experiments, the 
MACRO and Kamiokande results are not so well fitted. 
For the others, the agreement is excellent. 

The accelerator neutrino experiments are all well ac- 
counted for. It remains to be seen if further sampling of 
statistics in the LSND experiment will confirm its present 
value for P^ e . Finally, both reactor experiments are very 
well reproduced by the best point solution. 

The solution to the experimental constraints all give 
02 ~ 5° and 83 sa 45°. The mixing angle 8\ varies from 
15° to 45° as the flux factor r\ varies from 0.8 to 1.2 times 
its nominal value of 1. 

Since the CKM matrix can be written as U = 
^23(^1)^13(^2)^12(^3), where Uij(8) is a rotation by an 
angle 9 in the ij plane, it is clear that #3 = 45° corre- 
sponds to maximal mixing between the two lightest neu- 
trinos, v\ and V2- There are two particularly simple "so- 
lutions" more or less in the ranges of the mixing angles: 
(a) bimaximal mixing |23| with 0\ = 45°, 82 = 0, and 
83 = 45° and (b) single maximal mixing with 8\ = 30°, 
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62 = 0, and 63 = 45°. These two "solutions" both repro- 
duce data quite well (apart from the LSND experiment) 
the atmospheric ones with the appropriate factor r\ in- 
cluded. The CKM matrices for the simple "solutions" 
are given by 



U a = 




and 




to be compared with the CKM matrix for the best point 
solution 



U : 



0.7052 
-0.6441 
0.2964 



0.7052 0.0732 
0.5940 0.4820 
-0.3871 0.8731 



Note that the matrix element U e 3 is small compared to 
all the other matrix elements, in agreement with earlier 
suggestions p4| . For the simple "solutions," it becomes 
identically equal to zero. 

For the double ratio, we obtain approximately 



R = Rp/ Re 



1 RRfifx ~\~ Re.fi 



R Re, 



RP, 



(18) 



fie 



In our case, the CKM matrix is real so P efl = P^ e . Using 
the published |l^] values for the theoretical Vejv^ flux ra- 
tio for the upward-going muon-neutrinos (cosi? < —0.6), 
we have calculated the ratio of ratios for the solution 
obtained with 77 = 1. This ratio also eliminates the un- 
certainty in the muon-neutrino flux from the atmosphere. 
The result is 1Z — 0.51 ± 0.02, in good agreement with 
the value given by the Super-Kamiokande Collaboration 
of K = 0.41 ± 0.12 for multi-GeV data. This indicates 
that the electron-neutrino ratio R e is also in agreement 
with data. See Fig. || for the L/E dependence of 1Z. 

Our solution is also consistent with the (UP/DOWN)^ 
asymmetry measured by the Super-Kamiokande Col- 
laboration jl], as well as with the corresponding 
(UP/DOWN) e asymmetry. With our approximations, 
the (UP/DOWN)^ asymmetry is given by 



P, 



( jt 



P, 



nfi, 



(19) 



where is the number of /z-like upward-going events 
and is the number of /x-like downward-going events. 
The measured value for multi-GeV events is U^/D^ — 
0.54 ± 0.07 [Q, which is consistent with the value P w = 
0.56 ± 0.07 obtained from our analysis when r\ = 1. 

In the present analysis, we have deliberately considered 
mass ranges for the solution that avoid the region where 



the solar (and atmospheric) neutrinos could be affected 
by the MSW effect. Wc plan to address this question 
specifically in a future analysis. The mass range for the 
solar neutrino problem will then interfere with the anal- 
ysis of atmospheric neutrinos. 
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TABLE I. Neutrino experiments. The abbreviations SBL and LBL stand for short-baseline and long-baseline, respectively. 



Expcrirn. 


^ 


Type 


Reaction 


L (m) 


E (McV) 


E/L (oV 2 ) 


p 



GALLEX | 
Home-stake j 
Supcr-Kamiok 
Kami oka 



gdc 



• an d o 

lr 



solar 
solar 
solar 
solar 



1.496 X 10 
1.496 x 10 
1.496 x 10 
1.496 X 10 



n 



1.1 

i 5.5 
i 10 
J 10 



10" 
10" 



10" 



i2_ 10 -io 
n_ 10 -lo 

io-^-io- 10 

n-ll_ 10 -10 



0.60 ± 0.06 
0.33 ± 0.029 a 
0.474 ± 0.020 a 

0.54 ±0.07 ' 



Soudan 

Super-Kaminkandc 
MACRO ]3j,_, 
Kamiakandc 
1MB [ReL, 



atmospheric 
atmospheric 
atmospheric 
atmospheric 
atmospheric 



1.03 X 10' 
1.03 x 10 7 
1.22 x 10 7 
1.22 x 10 7 
1.03 x 10 7 



=s 1 000 
a 10 000 
i 100 000 
i 100 000 
i 250 000 



~ io-* 

~ 10~ 3 
~ 10~ 2 
~ 10 -2 
10 _1 -10 2 



~ 0.45 ±0. IE 
0.56 ± 0.07 
~ 0.5 ± 0.1 ]f5 

~ 0.85 ± 0.15 
~ 0.7 ± 0.1 [|17r 



LSNDTL9J 


1 — j accelerator (SBL) 




— ► V e 


30 


48 


1-2 


0.0026 ± 0.0010 £11.0005 


KARMBK 


r-fc'l accelerator (SBL) 


v» 


— * U & 


17.7 


a 40 


1-4 


< 0.0031 


NOMAD [ 


2J}T accelerator (SBL) 


v„. 


— > y r 


625 


24000 


~ 10 


< 0.0021 


CHORUS 


JEzfl accelerator (SBL) 


u t- 


— * Vt 


600 


26620 


10-100 


< 0.0006 


CHOOZ_( 


6| reactor (LBL) 


Ve 


— ► Ue 


1030 


3 


~ 10 -3 


0.98 ± 0.04 ± 0.04 


Bugey |7 


reactor (SBL) 




— > l> e 


40 


3 


10 _2 -10 _1 


0.99 ± 0.01 ± 0.05 



"Normalized according to the Bahcall-Basu-Pinsonneault 1998 standard solar model (BBP98 SSM) [pp^ . 



TABLE II. Probabilities. The best point solution obtained in our simulation, using first stochastic (s.) and then deterministic 
(d.) minimization, is compared to other "solutions." For comparison of the different solutions, the \ 2 values have also been 



listed. 



Experiment 




-Pexp 


P (best; s. + d.) 




P a 


P b 


SAGE 




0.52 ± 0.06 


0.49 




0.50 


0.50 


GALLEX 




0.60 ± 0.06 


0.49 




0.50 


0.50 


Homcstake 




0.33 ± 0.029 


0.49 




0.50 


0.50 


Supcr-Kamiokandc 




0.474 ± 0.020 


0.49 




0.50 


0.50 


Kamiokando 




0.54 ± 0.07 


0.49 




0.50 


0.50 


Soudan 




~ 0.45 ± 0.15 


0.45 




0.27 


0.35 


Supcr-Kamiokandc 




0.56 ± 0.07 


0.56 




0.45 


0.54 


MACRO 




~ 0.5 ± 0.1 


0.64 




0.50 


0.62 


Kamiokando 




~ 0.85 ± 0.15 


0.64 




0.50 


0.62 


1MB 




~ 0.7 ± 0.1 


0.64 




0.50 


0.62 


LSND 




0.0026 ± 0.0015 


0.0030 




3.8 x 10 _s 


4.2 x 10~ B 


KARMEN 




< 0.0031 


0.0017 




1.9 x 10~ 8 


2.1 x 10~ 8 


NOMAD 




< 0.0021 


9.5 x 10" 4 




1.1 x 10~ 3 


9.9 x 10~ 4 


CHORUS 




< 0.0006 


7.1 x 10~ 4 




8.2 x 10 -4 


7.4 x 10~ 4 


CHOOZ 




0.98 ± 0.08 


0.97 




0.98 


0.98 


Bugey 




0.99 ± 0.06 


0.99 




1.00 


1.00 


X 2 = F 






3.48 x 10~ 2 




6.22 x 10~ 2 


6.14 x 10~ 2 


a Parameter values: 


0i =45°, 


02 = 0, 3 = 45°, Am 2 


= 3.5 x 1(T 4 eV 2 , AM 2 


= 1 eV 2 , and r) = 1.2. 




b Parameter values: 


0i = 30°, 


2 = 0, 3 = 45°, Am 2 


= 3 x 1CT 4 eV 2 , AM 2 = 


1.1 eV 2 , and 77 = 1.0. 
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a) 

50 | , , , , , r 




2e-04 0.5 

0.7 0.8 0.9 1.0 1.1 1.2 1.3 0.7 0.8 0.9 1.0 1.1 1.2 1.3 

FIG. 1. The dependence on rj of the mixing angles (a) and mass squared differences (b) and (c). For all data points except 
rj — 1.0, we have used N = 10 instead of N = 20. 




1(T 2 10"' 10° 10 1 10 2 10 3 10 4 10 5 
I (eV 2 ) 



FIG. 2. The "ratio of ratios" 1Z as a function of / ~ 1.27L/E. The data points are obtained from Fig. 4 in Ref. Q. The 
dotted curve is a theoretical curve using the best point solution. The dashed and solid curves are obtained from the theoretical 
curve by using running average with different averaging lengths. The solid curve uses a larger averaging length than the dashed 
curve. 
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